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Vertical Structure of Neutrino-Dominated Accretion Disk and 
Applications to Gamma-Ray Bursts 

Tong Liu^, Wei-Min Gu^, Zi-Gao Dai^, and Ju-Fu Lu^ 
ABSTRACT 

We revisit the vertical structure of neutrino-dominated accretion flows in 
spherical coordinates. We stress that the flow should be geometrically thick 
when advection becomes dominant. In our calculation, the luminosity of neutrino 
annihilation is enhanced by one or two orders of magnitude. The empty funnel 
along the rotation axis can naturally explain the neutrino annihilable ejection. 
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Introduction 



> 

^sD ! Gamma-Ray Bursts (GRBs) are short-lived bursts of gamma-ray photons occurring at 

cosmological distances. GRBs are usually sorted of two classes (Kouveliotou et al. 1993): 
I short-hard GRBs (Tqq < 2s) and long-soft GRBs (Tgg > 2s). The likely progenitors are 
CN ■ the merger of two neutron stars or a neutron star and a black hole (Eichler et al. 1989; 

Q\ \ Paczyhski 1991; Narayan et al. 1992) and collapsar (Woosley 1993; Paczyhski 1998), re- 
spectively. The popular model of the central engine, namely neutrino dominated accretion 
flows (NDAFs), involves a hyperaccreting black hole with mass accretion rates in the range 
^ I of 0.01 ~ 10MqS~^. Such a model has been widely investigated in the past decade (see, 
; e.g., Popham et al. 1999; Narayan et al. 2001; Kohri & Mineshige 2002; Di Matteo et al. 
2002; Rosswog et al. 2003; Kohri et al. 2005; Lee et al. 2005; Gu et al. 2006; Chen & 
Beloborodov 2007; Liu et al. 2007; Kawanaka & Mineshige 2007; Janiuk et al. 2007; Liu et 
al. 2008). The model can provide a good understanding of both the energetics of GRBs and 
the processes of making the relativistic and baryon-poor fireballs by neutrino annihilation or 
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magnetohydrodynamic processes (see, e.g., Popham et al. [1999] and Di Matteo et al. [2002] 
for references). 

In cylindrical coordinates {R, z, (p), Gu & Lu (2007) discussed the potential importance 
of taking the explicit form of the gravitational potential for calculating slim disk (Abramowicz 
et al. 1988) solutions, and pointed out that the Hoshi form of the potential (Hoshi 1977), 

^^;{r,z)^^|J{r,0) + ^nlz\ (1) 

is valid only for geometrically thin disks with H/R < 0.2. Thus the well-known relationship 
Cs/^IrH — constant does not hold for slim disks with H/R < 1, where Cg is the sound 
speed, and Qk is the Keplerian angular velocity. Moreover, with the explicit form of the 
gravitational potential, Liu et al. (2008) found that NDAFs have both a maximal and a 
minimal possible mass accretion rate at their each radius, and presented a unified description 
of all the three known classes of optically thick accretion disks around black holes, namely 
Shakura-Sunyaev disks (Shakura & Sunyaev 1973), slim disks, and NDAFs. These works 
are, however, based on the following simple vertical hydrostatic equilibrium: 

i^ + ^ = o. (2) 

poz OZ 

instead of the general form (Abramowicz et al. 1997): 

p dz dz dR ^ dz 

where p is the mass density, p is the pressure, is the cylindrical radial velocity, and 
is the vertical velocity. Since Vz is not negligible for geometrically thick or slim disks, the 
solutions in Gu & Lu (2007) and Liu et al. (2008) are still not self-consistent. Recently, 
Gu et al. (2009) revisited the vertical structure in spherical coordinates and showed that 
advection-dominated accretion disks should be geometrically thick rather than being slim. 
However, the detailed radiative cooling was not considered in that work, and therefore no 
thermal equilibrium solution was established. 

The purpose of this paper is to investigate the vertical structure of NDAFs with detailed 
neutrino radiation. In section 2, with the self-similar assumption in the radial direction, we 
numerically solve the differential equations of NDAFs in the vertical direction. In section 3, 
we present the vertical distribution of physical quantities and show the geometrical thickness 
and the energy advection of the disk. In section 4, wc estimate the luminosity of neutrino 
annihilation and discuss some applications to GRBs. Conclusions are made in section 5. 
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2. Equation 

Wc consider a steady state axisymmetric accretion flow in spherical coordinates (r, 9, 
0), i.e., d/dt = d/d(f) = 0. We adopt the Newtonian potential ip = —GM/r since it is 
convenient for self-similar assumption, where M is the mass of the central black hole. The 
basic equations of continuity and momentum are the following (see, e.g., Xue & Wang 2005; 
Gu et al. 2009): 



1 d 1 d 



dvr , ve.dvr . GM Idp , . 

+ —{-KE- - ^e) = 5 (5) 

or r oo r p or 



+ —{^ +Vr) - —cote = — , 6 

Or r oO r pr oO 



H ^ + —{vr + vecote) ^ ——{r-'Tr4>), (7) 

r oO r pr-^ Or 



dr 



where Vr, ve, and are the three components of the velocity. Here, we only consider the 
r0-component of the viscous stress tensor, T^^ = pvrd{v^/r)/dr. The kinematic coefficient 
of viscosity takes the form: v — ckc^/Qk (e.g., Narayan & Yi 1995), where the sound speed 
Cs is defined as = p/p, the Keplerian angular velocity is Qk = {GM/r^Y^^, and a is a 
constant viscosity parameter. 

To avoid directly solve the above partial differential equations, some radial simplifica- 
tion is required since our main interest is the vertical distribution. Based on the radial self- 
similar assumption, Begelman & Meier (1982) studied the vertical structure of geometrically 
thick, optically thick, supercritical accretion disks. Under the same self-similar assump- 
tion, Narayan & Yi (1995) investigated the vertical structure of optically thin advection- 
dominated accretion flows (ADAFs). In fact, since the well-known self-similar solutions of 
ADAFs (Narayan & Yi 1994), such type of solutions has been widely investigated for different 
classes of accretion, such as slim disks (Wang & Zhou 1999), convection-dominated accre- 
tion flows (Narayan et al. 2000), NDAFs (Narayan et al. 2001), and accretion flows with 
ordered magentic fleld and outflows (Bu et al. 2009). Even though the detailed radiation 
was considered in some works (e.g., Di Matteo et al. 2002; Chen & Beloborodov 2007), and 



-4- 



therefore the solutions cannot be regarded as self-similar solutions, the self-similar assump- 
tion was still adopted such that the original differential energy equation can be simplified 
as an algebraic one. Furthermore, for optically thick flows, Ohsuga et al. (2005) showed 
that their simulations arc close to the self-similar solutions of the slim disk model (e.g., the 
density profile in their Fig. 11 ). In our opinion, the radial simplification is necessary for 
the study of vertical structure and it is a good choice to take the well-known self-similar 
assumption. 

Similar to Narayan & Yi (1995), we adopt the following radial self-similar assumption: 

p(r, e) oc r-^/^ (8) 



Cs{r,e),Vr{T,9),v^{r,e)<XT-^'\ (9) 

ve{r,e)^Q. (10) 
With the above assumption, equations (5-7) can be simplified as follows: 

+ ^c,2 + V - r V = 0, (11) 



-^ = Vcote, (12) 



Vr — (13 

2rf]K 

Integrating equation (4) over angle we obtain the mass accretion rate, 

M = -47rr^ / ' fWr sin OdO, (14) 

where 6*0 is the polar angle of the surface. 
The equation of state is 

P = Pgas + Prad + Pe + Pu, (15) 
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where Pgas, Prad, Pe, ^nd Pi, are the gas pressure from nucleons, the radiation pressure of 
photons, the degeneracy pressure of electrons, and the radiation pressure of neutrinos, re- 
spectively. Detailed expressions of the pressure components were given in Liu et al. (2007). 
We assume a polytropic relation in the vertical direction, p — Kp'^/^ , where is a constant. 

The energy equation is written as 

Qvis = Qadv + Qui (16) 

where Qvis; Qadv; ^^d Qi/ are the viscous heating rate per unit area, the advective cooling rate 
per unit area, and the cooling rate per unit area due to the neutrino radiation, respectively. 
Here we ignore the cooling of photodisintegration of a-particles and other heavier nuclei. The 
viscous heating rate per unit volume q^^^ = vpr'^[d{v^/r) / dr]"^ and the advective cooling rate 
per unit volume gadv = pVr{de/dr — p/ p^dp/dr) (e is the internal energy per unit volume) 
are expressed in the self-similar formalism as 

_ 3 (p-pe)^r , 
Qadv — —7; ■ [^°) 

2 r 

where the entropy of degenerate particles is negligible. Thus the vertical integration of Qvis 
and Qadv are the following: 

<5vis = 2 / ' qyi^rsmOde , (19) 



Qadv = 2 / Qsdyr sin 9d9. (20) 

The cooling due to the neutrino radiation Q^, can be written as 

= 2 / ' q^rsinOde , (21) 

where qi, is the sum of Urea processes, electron-positron pair annihilation, nucleon-nucleon 
bremsstrahlung, and Plasmon decay (see, e.g., Liu et al. 2007). We therefore can obtain the 
luminosity of neutrino radiation L^, by integrating Q,^. 

In our system, we have six physical quantities varying with 9, i.e., Vr, Vfj,, Cg, p, p, and 
T. The six equations for solving these quantities are equations (11-13), (15), the polytropic 
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relation, and the definition of Cs (c^ — p/ p). Three boundary conditions are required to solve 
the system since there is one differential equation, and the boundary 9q and the constant 
parameter K in the polytropic relation are unknown. Now we have already two boundary 
conditions, i.e., equations (14) and (16), thus one more boundary condition is required for 
solving the system, which is set to be = (accordingly p = and p = 0, e.g., Kato ct al. 
2008, p. 244) at the surface of the disk, i.e., 9 = 9q. The numerical method is as follows. For 
given a, M, M, r, and a test from the above six equations and two boundary conditions 
(except the energy equation, Eq. [16]), we can numerically obtain the vertical distribution 
of the above six quantities. With equations (19-21), we then check whether equation (16) is 
satisfied for the test ^o- By varying 9q we can find the exact value of for which equation 
(16) is matched, and therefore we obtain the exact vertical distribution of all the variables. 
In our calculations we take a = 0.1 and M — SM©. 



3. Numerical Results 

Figure 1 shows the variations of the density p, temperature T, electron fraction Y^, 
and radial velocity Vr with the polar angle ^ for M = 1MqS~^. Here is defined as 
= rip/ (rip + nn), where np and rin are the total number density of protons and of neutrons, 
respectively (e.g., Beloborodov 2003; Liu ct al. 2007). The solid, dashed, and dotted lines 
represent the solutions at r/vg = 10,40, and 100, respectively. The profiles of p and v,. 
are similar to that of the optically thin advection-dominated accretion flows (Narayan & Yi 
1995), i.e., p and Vr (the absolute value) decrease from the equatorial plane to the surface. 
On the contrary, electron fraction Y^ increases from the equatorial plane to the surface and 
approaches 0.5 near the surface, which means that the matter is non-degenerate. The vertical 
distribution of Vr, as shown in Fig. l{d), indicates a multilayer flow with the matter close to 
the equatorial plane being accreted much faster than that near the surface. 

Figure 2(a) shows the variation of the half-opening angle of the disk (7r/2 — ^o) with 
radius r/r^, where Vg — 2GM/c^ is the Schwarzschild radius. The solid, dashed, and dotted 
fines represent the solutions with M/MqS~^ — 0.1, 1, and 10, respectively. It is seen that, in 
the inner region of the disk, the half-opening angle increases as increasing accretion rates. 
For M = 1OM0S^^, the inner disk is extremely thick with the half-opening angle is ~ 1.4 
radian, which implies tfiat tfiere exists a narrow empty funnel ~ 20° along tfie rotation 
axis. Figure 2(6) sfiows tfie variation of tfic energy advcction factor /adv (= Qadv/Qvis) 
witfi r/vg. It is seen tfiat advcction becomes important in tfie inner disk for M > IMqS^^. 
Comparing Figs. 2(a) and 2(6), we find that the curves of the half-opening angle and the 
advcction factor are similar, which indicates that the geometrical thickness is relevant to 
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the advection. For f^dv — 0.5, it is seen from Fig. 2 that the half-opening angle is around 
1.3 radian. We therefore stress that NDAFs should be significantly thick when advection 
becomes dominant, which is in agreement with Narayan & Yi (1995) since their solutions 
imply that the flows are extremely thick with the half-opening angle approaching tt/2. 

4. Applications to GRBs 

In our calculations, the inner disk will be quite thick for large mass accretion rates, 
M > IMqS""^. Thus the volume above the disk shrinks and the radiated neutrino density 
increases. Accordingly, the neutrino annihilation efficiency also increases. We have obtained 
the neutrino luminosity L^, (before annihilation), thus the luminosity of neutrino annihilation 
L,yp can be roughly evaluated by the assumption: 1] oc V^~^ (sec, e.g., Mochkovitch et al 1993), 
where rj = Lyp/ Ly is the annihilation efficiency, and \4nn is the volume above the disk. For a 
given outer boundary rout) we calculate V^nn by integrating the region of ^ < a-nd r < Tout- 
The variations Lj, and L^p with M are shown in figure 3. The solid lines correspond to the 
present solutions whereas the dashed lines correspond to those in Liu et al. (2007). As shown 
in Fig. 3, for the same M, L^, is comparable, whereas L^y in the present results is significantly 
larger than that in Liu et al. (2007) by one or two orders of magnitude. Moreover, we find 
that for M = 5MqS-\ L is very close to L^^ which means that the density of radiated 
neutrino is so large that the annihilation efficiency is close to 1. Thus we can expect that 
Lyy is roughly equal to Ly for M > 5M0S~^. 

Many previous works have calculated the annihilation luminosity and claimed that the 
NDAF mode can provide enough energy for GRBs. However, GRBs are generally believed 
to be a jet with a small opening angle ^jct- The problem is that, the annihilation could not 
be limited into such a small angle even though the region well above the inner disk have 
larger luminosity than other place. We argue that our model is preferably to explain the 
ejection-like radiation, because the disk is adequately thick and there exists a narrow empty 
funnel along the rotation axis, which can naturally explain the neutrino annihilable ejection. 

5. Conclusions 

In this paper we revisit the vertical structure of NDAFs in spherical coordinates. The 
major points we wish to stress are as follows: 

1. We show the vertical structure of NDAFs and stress that the flow should be significantly 
thick when advection becomes dominant. 
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2. The luminosity of neutrino annihilation is enhanced by one or two orders of magnitude. 

3. The narrow empty funnel (~ 20°) along the rotation axis can naturally explain the 
neutrino annihilable ejection. 
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Fig. 1. — Variations of the density p, temperature T, electron fraction Y^, and radial velocity 
Vr with the polar angle 6, for which the given parameters are M/MqS~'^ = 1 and r/vg = 10 
(solid lines), 40 (dashed lines), 100 (dotted lines). 
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Fig. 2. — Variations of the half-opening angle of the disk (7r/2 — ^o) and the advection 
factor /adv with radius r/vg, for which the given parameter is M/MqS"^ — 0.1 (solid hne), 1 
(dashed hne), 10 (dotted hne). 



-13- 




Fig. 3. — Neutrino luminosity (thick lines) and annihilation luminosity L^p (thin lines) 
for varying mass accretion rates M . The solid lines correspond to the present solutions, 
whereas the dashed lines correspond to the solutions of Liu et al. (2007). 



